MAX-PLUS OBJECTS TO STUDY THE COMPLEXITY OF GRAPHS 



CRISTIANO BOCCI*, LUCA CHIANTINit, AND FABIO RAPALLO* 

Abstract. Given an undirected graph G, we define a new object Hq, called the mp-chart of G, in 
the max-plus algebra. We use it, together with the max-plus permanent, to describe the complexity 
of graphs. We show how to compute the mean and the variance of Hq in terms of the adjacency 
matrix of G and we give a central limit theorem for Hq . Finally, we show that the mp-chart is easily 
tractable also for the complement graph. 

' Key words, permanent of adjacency matrices; combinatorial central limit theorem; random 

(N ■ permutations; complement graph. 
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1. Introduction. The work presented in this paper has been inspired by the 
need of simple and actual techniques to measure the complexity of a graph, especially 
in the case of sparse graphs. This problem arises in several fields of applications, 
I from Computer Science to Economics, from Biology to Social Sciences. As general 

■ references for the graph theory, we mention in particular the books [7], [5], [S] and 
(-H I |10] . where the reader can find the main mathematical achievements in the theory. 

For a general survey on recent applications of graph theory see, for instance, [T]. The 
reader interested in some more technical papers can refer to |16| for applications to 
Economics, to [9] for applications to Econophysics, to [17] for applications to Biology, 
and to |15j for applications to Molecular Biology. In such papers, sparse graphs play 
PsJ ■ a prominent role. 

\ Although an undirected graph G = {V, E) is a rather simple structure, consisting 

of a set V oi N vertices and a set of edges E C V xV, \n graph theory there are 
several different approaches, depending on the specific application we are looking for. 
In particular, a graph can be fixed or random, depending on whether the elements in 
E are random or not. Moreover, many efforts have been done to analyze dynamical 
graphs, where the vertex set V and/or the edge set E vary with time, see e.g. [TT] . 

Here, we restrict our analysis to fixed graphs. In this framework, there are inter- 
esting developments in the area of Combinatorics, about the study of the properties of 
0—1 matrices. These matrices naturally arise in the framework of graphs, as the ad- 
jacency matrix Aq of the fixed graph G. Some recent developments in this direction, 
. , with applications to graph theory, are described in [5], [5], and [1]. 

■ In the present paper, we investigate some questions about undirected graph, in 
order to study and describe their structure, with special attention to sparse graphs. 
Our work is related to the matching problem. As a preliminary remark, we argue that, 
for sparse graphs, the classical descriptors of the complexity, such as the degree dis- 
tribution and the permanent of the adjacency matrix do not give actual information. 
Thus, we use the max-plus arithmetic and the corresponding expression of the perma- 
nent, and we show that this object is more suitable for describing of the complexity 
for sparse graph. The use of the max-plus arithmetic naturally leads to the definition 
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of a more complete index of the structure of a graph, and we define a vector cahed the 
mp- chart of the graph. This vector is nothing else but a probability distribution, and 
we show that it converges to a Normal distribution through a combinatorial Central 
Limit Theorem. Several examples on small- and medium-sized graphs are given to 
show that our definitions are easy to apply and provide practical information about 
the complexity of the structure of the graph under study. All the computations have 
been carried out with Maple, see [18], and R, see [19 . All the simulations come from 
simple R routines, without any additional package. 

This paper is only concerned with undirected fixed graphs. Nevertheless, the 
same strategy can be applied to other situations, such as bipartite graphs, undirected 
unfixed graphs, random graphs, and so on. 

The paper is organized as follows. In Section[2]we define the max-plus permanent 
of an adjacency matrix (i.e., the permanent under the max-plus arithmetic), we state 
its main properties, we study its connections with the classical permanent, and we 
discuss some simple examples. In Section |3] we define a new object associated to a 
graph, and we call it the mp-chart of the graph. We compute its mean and variance, 
and we show that, under suitable conditions, it converges to a Normal distribution 
as the size of the graph goes to infinity. Some simple simulations show that the 
convergence is quite good also for small values of the size. In Section HI we show how 
the mp-chart of a graph is related to the mp-chart of the complement graph. Finally, 
Section [5] is devoted to suggest some future directions of this research. 

2. The max-plus permanent. Let G = {V,E) be an undirected graph with 
N vertices. Let Ac be the N x N adjacency matrix of G, defined by (^G)ij = 1 if 
e E and otherwise. In the classic definition of undirected graph, the matrix 
Aq is symmetric and with zero diagonal entries, as we do not consider loops. 

As mentioned in the Introduction, the study of the complexity of a given graph 
is one of the most relevant problems about graphs in Applied Probability. This 
analysis can be performed through the distribution of the degrees (i.e., the number of 
edges involving each vertex) and through the permanent (or the determinant) of the 
adjacency matrix Aq. 

The determinant of Aq is 

N 

det(AG) = 5](-l)l'^in(^G).,.(.), 

TT i—1 

where the sum is taken over all the permutations tt of {1, ... , N} and |7r| denotes the 
parity of tt. The permanent of Aq is 

N 

7r i—1 

The use of permanent to describe the complexity of a graph is justified by the 
following well-known property. 

Proposition 2.1. The permanent of Aq is the number of bijections (j) : V ^ V 
compatible with E, i.e. such that (v,(j){v)) G E for all v G V. 

In fact, perm(yl(5) is the number of permutations tt with 7r(i) = . . . = A7v.7r(Af) = 
1 and the permutation tt is just the bijection (f) in the proposition. 

However, the analysis based on the degree distribution and the permanent is not 
adequate for sparse graphs. In fact, it is enough to have an isolated vertex to produce 
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a null permanent. Nevertheless, it is interesting to study the structure of a sparse 
graph. 

To overcome this difSculty, we make use of the tropicalization of the permanent. 
In the classical settings, Tropical Arithmetic is defined through the operations: 

X ® y — min{a;, y} x ® y — x ^- y 

But, with Tropical Arithmetic, the determinant (or permanent) of an adjacency 
matrix is always 0, because of the nullity of the main diagonal of Aq- 
Thus, we use the max-plus algebra, with operations: 

X ® y — max{a;, y} x ® y ~ x + y 

Consequently, the explicit expression of the max-plus permanent is 

CN \ N 

(^(^G)j,7r(j) j = max^(AG),,^(j) (2.1) 

The max-plus permanent is the maximum over N\ terms. Each of the N\ terms 
is the sum of N terms in {0, 1}. Thus, the max-plus permanent permi„p(^G) is zero 
if and only if the matrix Aq is the null matrix. On the opposite side, the maximum 
allowed value of the max-plus permanent is TV. 

The use of the max-plus permanent to analyze sparse graphs has a first reason in 
the following property. 

Lemma 2.2. The following relation holds: 

permmp(AG) ^ perm(AG) > . (2.2) 

Proof. permmp(AG) = A'^ if and only if there exists a permutation tt such that 
'Y^^=i{AG)i^T![i) — N . This happens if and only if there exists tt such that (Ag)^ = 
1 for i = 1, . . . , A^, i.e. if and only if perm(AG) > 0. □ 

Remark 2.3. Notice that, from Lemma \ and from the previous discussion, 
it follows that the max-plus permanent is able to discriminate among graphs with 
standard permanent equal to zero. 

Moreover, we explicitly write the following consistency property, whose simple 
proof is straightforward. 

Lemma 2.4. Let G and H be two graphs on two disjoint sets of vertices. Then, 

permmp(AGuH) = permnip(AG) + perm„ip(A_f/) . (2.3) 

The max-plus permanent has interesting connections with the subgraphs. Let 
G" = {V',E') a graph. UV CV and E' C E, then G' is a subgraph of G = iV,E). 
A subgraph G' = {V, E') is the subgraph induced by V if E' contains all the edges in 
E involving the vertices in V . In order to analyze the max-plus permanent, in view 
of Equation (12. ip . we introduce here the notion of i-tcrm, which is strictly related to 
the subgraphs of G. Such connections will be studied later in this section. 

Definition 2.5. Given a graph G with adjacency matrix Ac, a t-term is a 
sequence of indices {ii,ji) ■ ■ ■ (*t, Jt) with 
• 1 < ii < . . . < it < N; 
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• the jk 's, with 1 < jj. < N are all distinct; 

• (^G)ifcjfc = 1 for all k. 

For a t-term P, we denote I{P) = {ii, . . . ,it} and J{P) = {ji, . . . , jt}. 

Roughly speaking, a i-term corresponds to a sequence of positions of t ones in 
the permutations. A straightforward consequence is the following statement. 

Proposition 2.6. The max-plus permanent of G is q if and only if there exists 
a q-term and there are no t-term with t > q. 

Proposition 2.7. Let t the maximum integer such that there exists a t-term, 
then there exists a t-term P such that I{P) — J{P)- 

Proof. We prove the statement by induction on t. If t = 2 there is nothing to 
prove since if the 2-term is given by (ii, Ji)(i2ii2) it is enough to consider the 2-term 

Consider now a i-term Q and suppose there exists a k such that i^ is in I{Q)\J{Q)- 
First of all we notice that jk must be in I{Q). If not, we can add the element {jk , ik) to 
Q obtaining a + l)-term which is a contradiction, since Q is maximal. Hence, since 
jk G I{Q) then there exists a s such that jk = is & I{Q)- Then, we substitute {is,js) 
with {jk,ik) in our i-term and we obtain a new t-term of the form {ik, jk){jk,ik)Q' 
where Q' is a (t — 2)-term. This term Q' arises from the sub-matrix A' of A{G) where 
we remove rows and columns ik and jk- Hence Q' is a maximal (t — 2)-term for A' . 
By induction, the proof follows. □ 

Remark 2.8. If P is a 3-term, then we must have I{P) = J{P)- In fact, if 
P is (ii, ji)(«2, J2)(*3, is) with I{P) 7^ JiP) then, by the previous proposition, we 
obtain a new 3— term {iki jk){jk,ik){i3T j^,)- Then it would be possible to extend it to 
{'ik,jk),{jk,ik){'i3,j3){j3,h) against the maximality of P. 

Remark 2.9. In view of Proposition the max-plus permanent is the cardi- 
nality of the largest subset of V with a bijection compatible with E. This is another 
way to see that the max-plus permanent is able to detect the complexity of the graphs 
with null classical permanent. 

Denote by £g the number of edges of a graph G. Among the subgraphs of Propo- 
sition 12. 7[ we are mainly interested in the ones with a minimal number of edges. 
These subgraphs are maximal in term of permmp ( Ac ) j but minimal in term of ^G'■ 
We made this more precise by the following definition. 

Definition 2.10. An mp-maximal subgraph G' of a graph G, is a subgraph of G 
with q = perminp(AG) vertices, 



permmp(AG') = permmp(AG) 



(2.4) 



and for all other subgraph G" of G satisfying \2.4^ one has ic < (g" ■ 

The rest of this section is devoted to the discussion of some examples and some 

useful remarks. In order to understand the definitions introduced above, we start 

with some small graphs. 

Example 2.11. Let us analyze the three graphs on 4 vertices drawn in Figure 

\2.1l Their adjacency matrices are respectively 



Ac 



/O 1 1\ 
10 10 
10 1 

yi 1 0/ 



Ag2 = 



/O 1 0\ 

10 

1 

\o 1 oy 



Ag. = 



/O 1 1 0\ 
10 10 
110 

\o oy 



In the first graph, all vertices are connected and permmp(AG'j) = 4. However, 
this is not the minimal way to obtain a max-plus permanent equal to 4. In fact, it is 
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Gi G2 G3 

Fig. 2.1. The graphs Gi, G2 artd G3 for Examvle \2.11\ 




Fig. 2.2. The two graphs Gi and G2 for Examvle \2.12\ 



easy to check that perniinp(AG'2) — 4. Thus the graph G2 represents a mp-maximal 
subgraph for the graph Gi, but it is not the only one. If we look now at the graph G3, 
we notice that there is a cycle of length 3 and an isolated vertex. In such case, we 
have pernimp(AG3) = 3, and there is only one mp-maximal subgraph. 

Example 2.12. To illustrate the behavior of the max-plus permanent and of the 
mp-maximal subgraphs, we analyze two opposite examples, with the same length. The 
two graphs are drawn in Figure \2. 2\ The graph Gi on the left is the union of a tree 
and two isolated vertices, with pernimp(AGj) = 2 and 4 maximal subgraphs with two 
vertices and one edge each. On the opposite side, the graph G2 has a perfect matching, 
perm„ip(AG2 ) = 6 and there is only 1 maximal subgraph, i.e., the graph G2 itself. 

Proposition 2.13. Two mp-maximal subgraphs are not disjoint. 

Proof Consider a graph G such that permmp(AG) — t and let G' — {V , E') and 
G" = {V" ,E") be two mp-maximal subgraphs of G with t vertices each. Suppose 
that V' and V" are disjoint. Then, by Formula ()2.3|) . the adjacency matrix of G'UG" 
has max-plus permanent 2t. Then permnip(AG) > 2t which is a contradiction. □ 

Remark 2.14. In the max-plus arithmetic, the definition of determinant is not 
unique, see JJ^. Therefore, one has to define the positive and negative determi- 
nant. In particular, the positive max-plus determinant is the maximum of the sums 

^i,7r(i) over all even permutations n. The negative determinant is defined by taking 
the odd permutations instead of the even ones. This issue is another reason to use 
the permanent instead of the determinant in the max-plus environment. 

3. The mp-chart of a graph. The information about a graph is not con- 
tained only in the max-plus permanent, but in the whole distribution of the iV! terms 
Thus, in this section we define the mp-chart of a graph as the dis- 
tribution of the A'^! terms above, and we prove that this distribution converges to a 
Gaussian distribution through the Hoeffding's combinatorial central limit theorem. 
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see [2]. We also show that the mean and the variance of that distribution can be 
computed easily from the adjacency matrix. 

Definition 3.1. Let G be a graph and Aq its adjacency matrix. Let hoik) he 
the number of permutations tt such that X]i=i(^G)i,7r(i) — We call the {N + 1)- 
dimensional integer vector Hq — {hQ{0), . . . , hclN)) the mp-chart of the graph G. 

This object captures many features of the graph and has some relevant theoretical 
properties. To understand the meaning of Hq^ notice that hcik) is just the number of 
permutations tt such that the sequence {Ag)i,tv{i), • • • , (^G)Af,7r(JV) contains a fc-term 
but not a (fc + l)-term. This gives precisely the meaning and the usefulness of the 
notion of random permutation in that context. 

Example 3.2. We use here a very simple scheme inspired by Econophysics, see 
\21}/ and 1 12^ . Consider a population with N agents, each possessing one good. The 
goods can be sent and received only along the edges of a graph G and each agent can 
possess only one good. Given a random permutation vr of {1, . . . ,N}, the i-th agent 
can send its good to 7r(i) if it receive a good from TT^^{i). The quantity X]iLi(^G)i,Tr(i) 
is exactly the number of agents involved in this process. Of course, similar examples 
can be adapted to many other sciences. 

We start the analysis of the mp-chart with the study of the mean E{Hg) and 
the variance Y{Hg). Although these computations could be carried out applying 
Theorem 2 in [14 , it is useful to state explicitly the proof for adjacency matrices. 

Given an adjacency matrix A — (Ag), note that, for a random permutation tt, 
the object 

5'A'(7r) = Ai^7r(l) + • ■ • + ^JV,7r(Ar) (3.1) 

is the sum of N binary random variables, and Tr(i) is the addendum chosen in the 
i-th row of the adjacency matrix. 

In order to analyze the mean of Hg, defined by 

mo) = ^^^^, (3.2) 

it is convenient to adopt an inductive approach. 

Theorem 3.3. The mean of the mp-chart Hg is 

E{Hg)^^, (3.3) 

where £g is the number of edges in the graph G. 

Proof. Clearly, if Ag is the null matrix, then K{Hg) = 0. Suppose that the 
formula p.3p holds true for — 1. By direct inspection, adding one edge has the 
following consequences. Among the iV! terms Sn{t^)- 

• {N — 2)1 of them increase by 2; 

• 2{N -2){N -2)1 oi them increase by 1; 

• the remaining {{N — 2)^ + N — 1)(-/V — 2)! do not change. 

Thus, 

^ 2{iG - 1) 2(jV-2)! + 2(jV-2)(jV-2)! ^ 2{Ig - 1) 2^2^ 
^ ^' N N\ N N N ' 

□ 

Example 3.4. Given a complete graph G, its adjacency matrix has on the 
diagonal and 1 elsewhere. The graph has N{N — l)/2 edges. Hence K{Hg) = iV — 1 
which is the maximum allowed. 
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Case (a) 



Case (5) 



Fig. 3.1. The two cases arising in the proof of Theorem 1 3. 51 



Notice that the mean E,{Hc) depends only in the number of edges of G, whatever 
they are collocated, that is, E,{Hg) does not take into account the topology of the 
graph. On the other hand, the variance V(i?G) depends on the position of the edges. 

Theorem 3.5. The variance of the mp-chart Hq is 



where di, . . . ^ (In are the degrees of the vertices and Tij =< ri,rj > is the scalar 
product of the i-th and the j-th row of Ac- 

Proof. By direct computation, the formula (|3.4|) holds for N < 2. 

To prove the validity of Eq. p.4|) for > 3, it is enough to compute the covari- 
ances 



Without loss of generality we can fix z = 1 and j ~ 2 and we write for brevity Ai for 
^i,7r(i) E^nd A2 for A2,7r(2)- Morcovcr, we suppose that di and c?2 are both non zero. 
(If'di = or da = 0, then trivially Cov(yli, A2) = 0). 

We divide the computation in two cases, and to help the reader we have sketched 
the two cases in Figure l3Tl 

• Case (a): (1, 2) is not an edge of the graph. Then: 



- there are 2{N — 1)! permutations such that 7r(l) = 1 or 7r(l) — 2. For 
all these cases, {Ai — 1,A2 = 1) is impossible; 

- there are 2{N-2){N-2)\ permutations such that 7r(l) > 2, but n{2) = 1 
or 7r(2) = 2. Also in all these cases, {Ai — 1,A2 = 1) is impossible; 

- there are {N - 2){N - 3){N - 2)! permutations such that 7r(l) > 2 and 
7r(2) > 2, and among these permutations 




(3.4) 



Cov(Aj_^(j),v4j,^(j)) = E{Ai^^(i)Aj^^(^j)) - didj/N^ = 

= P(^,^,(,) = 1, Aj- = 1) - ddj/N^ . 



{Ti,2d2 + (di - Ti,2)d2){N - 2)! = (did2 - Ti.2)(iV - 2)! 



are such that A1A2 — 1. 



Therefore, 



Cov(^i, A2) 



did2 - iVTi,: 
N^{N - 1) 
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Case (6): (1, 2) is an edge of the graph. Then: 

— there are {N — 1)1 permutations such that 7r(l) = 1. In ah such cases, 
(^1 = 1,A2 = 1) is impossible; 

— there are {N — 2)! permutations such that 7r(l) = 2 and 7r(2) = 1. For 
such permutations, Ai — 1 and A2 = 1] 

— there are (A*" — 2)(iV — 2)! permutations such that 7r(l) = 2 and 7r(2) > 2. 
Among these permutations, {d2 — 1){N — 2)! are such that A1A2 = 1. 

— there are (iV — 2)(A^ — 2)! permutations such that 7r(l) > 2 and 7r(2) = 1. 
Among these permutations, {di — 1){N — 2)! are such that A1A2 — 1. 

— there are (A^ — 2)(A^ — 2)! permutations such that 7r(l) > 2 and 7r(2) = 2. 
In all such cases, A1A2 = 0. 

— there are {N - 2){N - 3)(iV - 2)! permutations such that 7r(l) > 2 and 
7r(2) > 2, and among these 

(Ti,2(d2 - 2) + (di - 1 - Ti,2)(d2 - l))iN - 2)1 = 
= {did2 -di-d2~ ri,2 + l)(iV - 2)! 

are such that A1A2 — 1. 
Therefore, adding up all the contributions, we obtain again 



Cov(Ai, A2) 



did2 - NTi^2 

N^{N - 1) 



The formula in Eq. (13.41) is now straightforward. □ 
Example 3.6. Consider the matrices 



/o 


1 












/o 


1 










1 





1 










1 

















1 






















1 




























1 








Vo 











0^ 




\o 











0/ 



The graph Gi has two consecutive edges, while the graph G2 has two disjoint 
edges. An easy computation gives 



Hg, = (48,48,24,0,0,0) 



Hg2 = (53,44,18,4,1,0) 

with equal means E{Hgi) — E{Hg2) — 4/5. On the contrary, the variances are 
V(i/Gi) = 14/25 andYiHc^) = 19/25, respectively. 

Example 3.7. As a second example, consider the two graphs on the set vertices 
V — {1, . . . , 7} shown in Figure [KM Notice that Gi and G2 differ by only one edge. 
The two mp-charts Hgi and Hg2 have the same mean and variance, namely 

12 170 

E{Hg, ) = ) = y V(i?G, ) = V(i?G. ) - 

but the mp-charts are different: 



Hg, = (678, 1512, 1716, 840, 294, 0, 0, 0) 
Hg2 = (674, 1480, 1792, 840, 218, 32, 4, 0) . 
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The results above lead to a central limit theorem. 

Theorem 3.8. Let Gn be a graph with N vertices and let Aq,^ be its adjacency 
matrix. Let tt be a random permutation of {1, ... , N} chosen with uniform probability 
and define 



N 



Sn{tt) ^^{AGt,)i.yr(i) ■ (3.5) 



i=l 



//V(S'jv) goes to infinity as N oo, then the distribution of Sn is asymptotically 
normal. 

Proof. We make use of Theorem 3 in [14) . Define the auxiliary matrix R with 
elements 

R., - {Ac^k, - I - I + ]^ E(^Gj..fc (3.6) 

h,k 

Then, a sufficient condition for the asymptotic normality is that 

lim ^^^<^<^^. ^Q^ (3.7) 

AT Z^iJ=l ^'■i.j 

Now observe that the numerator is bounded, as — 2 < Ri^ < 2 for all i and j. 
Moreover, Theorem 2 in the same paper [T3] states that 

1 ^ 

Combining these facts, the result follows. □ 

Remark 3.9. Note that in our problem one can not use the classical central limit 
theorems based on a-mixing sequences or m-dependent variables, see for instance 
Ch. 27] and i20\j . Indeed, the covariance between ^i^7r(i) o,i^d ^j.7r(j") does not vanish 
''■s \i ~ j\ goes to infinity. 

In order to inspect the behavior of the convergence to the Gaussian distribution, 
we have computed the mp-chart for some graphs with 20 vertices. 
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Chart Distribution function 




5 10 15 20 2 4 6 8 10 12 



Fig. 3.3. Checking the convergence for N = 20; First graph and its results. 

The three examples in Figures I3.3II3.5I show that the convergence is quite good, 
meaning that the Gaussian approximation is vahd also for medium-sized graphs. 

Two remarks are needed to understand the examples: (a) The mp-chart is ap- 
proximated through a standard Monte Carlo technique, sampling 100,000 random 
permutations. This number is considerably smaller than the total number of permu- 
tations (20! = 10"'^^), but it provides quite accurate approximations; (6) The results 
are presented through two plots, showing the mp-chart (normalized to 1) and its 
distribution function, both compared with the appropriate Normal distribution. 

The first graph corresponds to an adjacency matrix with block structure. The 
graph and the two plots of the results are presented in Figure 13.31 

The second graph has a different shape, as it corresponds to an adjacency matrix 
with band structure. The graph and the two plots of the results are presented in 
Figure EH 

For the third graph analyzed here, we present only the results. The graph has 
been constructed with N{N — l)/4 = 95 edges randomly chosen among the 190 edges 
of the complete graph with uniform probability. The results are shown in Figure 13.51 

4. The mp-chart of the complement graph. In the literature, the comple- 
ment of a graph G = {V, E) is a graph on the same vertex set V and the set of edges 
V"^ \ E. Since our starting graph has no loop (e,e), this forces to contain all of 
them. To avoid this problem we give a different definition of complement graph, more 
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Fig. 3.5. Checking the convergence for N = 20: The results for the third graph. 
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useful for our purposes. 

Definition 4.1. Given a graph G = {V,E), its complement graph G'^ is a graph 
on the same vertex set V and the set of edges (V"^ \ /S) \ E, where A C E"^ is the 
diagonal set, i.e. A = : v G V}. 

Remark 4.2. From the previous definition, we notice that, if G and G"^ have 
respectively Ig and Ic- edges, then £g + ic = N{N — l)/2. 

As mentioned in Section [31 there are some nice properties linking the mp-chart 
Hg of a graph G with the mp-chart of its complement G"^. To study these connections, 
we start with a preliminary lemma. 

Lemma 4.3. Let G be a graph with N vertices. Denote by T[j the scalar product 
between the i-th row and j-th column in the adjacency matrix of G'^ . The following 
formula relates the quantities T^- and Tij : 

^ ^T,j+N-2-d,- dj+ < r„Ej > + <rj,E,> . (4.1) 

where the Ei 's are the vectors in the canonical basis of M.^ . 
Proof. Define the vector 

Vij ^ ri - r-j- < ri - r.j,Ei > E,- < ri - r.j,Ej > Ej (4.2) 

Since < ri — rj,Ei >=< ri,Ei > — < rj,Ei > and < rt,Et >= (in fact, the i-th 
coordinate of rt is zero, while Et has a 1 in the i-th coordinate and elsewhere), we 
can write 

Vij = Ti - r.j+ < r.j,Ei > E,~ < r,,Ej > Ej 

The scalar product < Vij, Vij > measures the number of positions, out of the diagonal, 
where ri and rj are different. Thus, if we denote by T[j the scalar product of the 
corresponding lines, , r| in the complement graph, one has 

nj - < v,j (4.3) 

Now, we substitute in the previous formula the expression of Vij given in (|4.2|) . and 
we obtain: 

T^j^ N -2-T,j- <v,j,v,j >=N-2-T,j- 

- <ri- rj+ < rj,Ei > E,~ < ri,Ej > Ej,ri - rj+ < rj,Ei > E,- < r,,Ej > Ej > . 

Noting that < rt,rt >= dt and < ri,rj >—< rj,ri >— Tij, a straightforward compu- 
tation leads to 

T,'j = + TV - 2 - - dj+ < r„Ej >2 + < rj,E, >^ . 

Since, for all i,j, with i ^ j, the value of < ri,Ej > can be either or 1, we can 
remove the squares from the previous formula, leading to Equation ()4.ip . □ 
Theorem 4.4. Given a graph G with N vertices, 

(a) E{Hg^)^N-1-E{Hg); 

(b) WiHG^) = YiHG) + 1 - 2EiHG)/{N - 1). 

Proof. To prove part (a) , it is enough to use Theorem 13.31 and Remark 14.21 One 

has 



MAX-PLUS OBJECTS FOR GRAPHS 



13 



To prove part (6), we apply Theorem [23] to the graph G"^. Therefore we have: 



i=l i,j=l 

The degree of a vertex in the complement graph is given hy d'^ — N ~ I ~ di. Using 
also Lemma 1113] one can write V(i/G'=) as 

i=l 

(N -l~d,)iN -1-dj)- NjTjj +N-2-d,-dj+<n,Ej> + < rj.Ej >) 
^ ^ 7V2(iV- 1) 

(4.4) 

The first sum in formula (14.41) can be written as 



A {d,){N -d,) A {N-\-2di) 

1=1 i=l 



where an easy computation shows that 

A (iV - 1 - 2d,) _ TV - 1 4£g 

About the second sum in formula (|4.4p . we observe that each term can be ex- 
pressed as 

(iV - 1 - d,)(jV - 1 - dj) - N{Tij +N-2-d^- dj+ < n,Ej > + <rj,Ei >) _ 

N^{N -1) ^ 
_ d,dj - NT,j +d, + dj + l- N{< ri,E^ > + < rj.Ei >) 
~ 7V2(7V_ 1) ■ 

Since 

Af TV 

^ < r„ i;, > + < r,, >) ^ Ug and ^ (d, + d,) = 4^G(iV - 1) 

the second sum becomes 

^ didj - NT.^j 4(iV - l)iG - N{N - 1) - 

N^{N-1) ^ 7V2(iV- 1) ■ 

Thus, we obtain 

4^G 



V(i/GO = V(i/G) + 1 



N{N - 1) 
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and, considering Theorem 13.31 the formula in (b) fohows. □ 

Some major remarks on the Theorem above are now in order. 
Remark 4.5. As a first trivial example, we consider the limit situation of an 
empty graph. Let G be the empty graph. Its mp-chart has K^Hg) = Y{Hq) — 0. 
In this case G"^ is the complete graph with N{N — l)/2 edges and, by Theorem \4-4\ 
part (b), one has W{IIqc) = 1. This can be verified also by direct computation. As 
a matter of fact, the adjacency matrix of G"^ consists of non-zero entries out of the 
diagonal. Thus d1 ~ N — \ for all i and T[j — N — 2 for all i,j , with i ^ j . Hence 



N 



i=l 



d^{N-d,) 

iV2 



djd^ — NTj 



N 

Yl^ N^(N -I) 



N-1 (TV - 1)2 - iV(iV - 2) 



N 



N^{N - 1) 



These computations show that the variance of the mp-chart of a complete graph is 
invariant on the number of vertices. 

Remark 4.6. Few straightforward algebraic calculations show that the difference 
W^Hg") ~Y{Hg) lies between —1 and 1. Therefore, under the hypotheses of Theorem 
\3.8\ when N goes to infinity we have thatY{HG<') = W{IIg). Intuitively, the difference 
between V(Hgc) and Y{Hg) depends on the diagonal entries which are forced to be 
zero in the adjacency matrix Ag" • The effect of these entries vanishes when the size 
of the graph goes to infinity. 

Remark 4.7. Another interesting property follows from Theorem \4-4\ First, 
notice that E(JfGc) = ¥.{Hg) implies that G and G'^ have the same number of edges, 
namely N(N — l)/4. (This is not possible for all values of N). In such a case, Hg" 
and Hg are forced to have the same variance, no matter how is complicated the graph 
G. 

The computation of the whole mp-chart of the complement graph G"^ from the 
mp-graph of G is less easy. Given a graph G, we build a {N + 1) x (N + 1)— matrix 
Mg, indexed, both on rows and columns, by {0, . . . ,N}, an defined as follows. The 
entry {MG)i,j is the number of permutations tt such that X]^i(^G)s,7r(s) — j and tt 
has i diagonal elements (that is, 7r(s) = s for i elements). 

Example 4.8. Consider the graph G with matrix 



Ag = 



^0 1^ 

10 

10 1 

yi 1 oy 



The corresponding matrix Mg is 



Mg = 



/I 2 

4 
3 



2 1\ 






\1 0/ 



The matrix Mg allows the computation of both the mp-charts Hg and Hg". 
Roughly speaking, to compute the mp-chart of G is is enough to sum the columns 
of Mg, while to compute the mp-chart of G"^ we need to sum the entries of suitable 
diagonals of Mg. More precisely, the following relations hold true. 

Proposition 4.9. For a graph G, we have for all j ~ 1, . . . , N : 
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(a) The components of the mp-chart Hq are: 



N 



(b) The components of the mp-chart Hqc 



are: 



ha- (j) 



N-j 



N- 



1=0 



The 



Proof. The first relation follows by the definition of mp-chart, as the sum of entries 
in the j-th column of Mq is the number of permutations with j elements equal to 1, 
that is hc{i)- 

To prove the second relation, it is enough to prove that for all i and j we have: 
{MGa ji j = {MG)i.N-j-i- Suppose that tt is such that X]^i(^G)s,7r(s) = j and tt has 
i diagonal elements. When we consider tt on A{G'^), we have A(G"^)s^^(s) = 1 for the s 
such that {Ag)s.tt{s) = 0, except for the diagonal entries where we still have 0. Hence 
there are N ~ j — i entries in ^(G"^) such that A{G'^)s,tt{s} = 1- This completes the 
proof. □ 

Example 4.10. Consider the complement graph G"^ of the graph G in Example 



/O 1 1 0\ 

10 1 

10 

\0 1 0/ 



/l 2 3 2 l\ 
4 4 
Mgc = 3 3 

\1 0/ 

We notice that in this simple case Mq = Mq^. This is due to the fact that G and G'^ , 
up to the labels of the vertices, are equivalent. 

Remark 4.11. In principle, the max-plus permanent of the complement graph 
G'^ can be computed from the matrix Mq- In fact. Proposition \4.^ shows that the 
mp-chart of the complement graph can be computed from the matrix Mq by adding 
along suitable diagonals and the max-plus permanent is just the position of the last 
non-zero element of the mp-chart. However, as the matrix Mq is not easy to compute 
for large graphs, this approach does not help in actual computations. 

5. Future directions. The max-plus permanent and the mp-chart studied in 
this paper lead to several new questions. In fact, we have analyzed here only the 
case of undirected graph. Therefore, among the future directions of our research, 
there will be the extension of the definition of max-plus permanent and mp-chart for 
undirected and weighted graphs. Moreover, special classes of graphs may be studied, 
such as bipartite graphs, or fixed-degree graphs. In particular, fixed-degree graphs 
correspond to adjacency matrices with fixed margins and in that context algebraic 
and combinatorial methods have demonstrated already their potential. 
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Another possible research direction is strictly in graph theory. As a matter of 
fact it could be interesting to compare the mp-chart of a graph with other well-known 
descriptors of its complexity. In a recent work in progress the mp-chart is compared 
to matching polynomials. By several examples we know that there exist different 
graphs with the same mp-chart. We notice that, in all these cases, also the matching 
polynomials coincide. So a principal question would be: If two non-isomorphic graphs 
have the same mp-chart, then their matching polynomials are equal? 

Finally, applications to large graphs, possibly through simulation techniques, will 
be investigated in order to use the tools presented in this paper to real data examples. 
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